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The Hi 21 cm absorption optical depth and the N(Hi) column density derived from Lya 
-i_ > ' absorption can be combined to yield the spin temperature (T s ) of intervening damped Lya ab- 

sorbers (DLAs). Although spin temperature measurements exist for samples of DLAs with 
Zsbs < 0.6 and z a b s > 1-7, the intermediate redshift regime currently contains only two 
Hi 21 cm detections, leading to a "redshift desert" that spans 4 gigayears of cosmic time. 
To connect the low and high redshift regimes, we present observations of the Lya line of six 
0.6 < z a bs < 1.7 Hi 21 cm absorbers. The dataset is complemented by both low-frequency 
Very Long Baseline Array (VLBA) observations (to derive the absorber covering factor, /), 
and optical echelle spectra from which metal abundances are determined. Our dataset there- 
fore not only offers the largest statistical study of H 1 21 cm absorbers to date, and bridges the 
redshift desert, but is also the first to use a fully /-corrected dataset to look for metallicity- 
*0 ■ based trends. 

The metallicities of Hi 21 cm absorbers are consistent with those of optically-selected 
DLAs at the same redshift. In agreement with trends found in Galactic sightlines, we find that 
the lowest column density absorbers tend to be dominated by warm gas. In the DLA regime, 
spin temperatures show a wider range of values than Galactic data, as may be expected in a 
heterogenous galactic population. However, we find that low metallicity DLAs are dominated 
by small cold gas fractions and only absorbers with relatively high metallicities exhibit sig- 
| nificant fractions of cold gas. Using a compilation of H i 21 cm absorbers which are selected 

to have covering-factor corrected spin temperatures, we confirm an anti-correlation between 
metallicity and T, at 3.4<x significance. Finally, one of the DLAs in our sample is a newly- 
discovered Hi 21 cm absorber (at z a b s = 0.602 towards J 143 1+3952), which we find to have 
the lowest /-corrected spin temperature yet reported in the literature : T s = 90 + 23 K. The 
observed distribution of T, and metallicities in DLAs and the implications for understanding 
the characteristics of the interstellar medium in high redshift galaxies are discussed. 

Key words: galaxies: high redshift, galaxies: ISM, galaxies: abundances, quasars: absorption 
lines. 



* Based on observations made with ESO Telescopes at the La Silla Paranal 
Observatory under programme IDs 082.A-0569, 166.A-0106, and 085.A- 
0258, and observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained from the data archive at the Space Telescope Institute. 
STScI is operated by the association of Universities for Research in As- 



1 INTRODUCTION 

The interstellar medium (ISM) gas in galaxies represents the raw 
material for star formation through cosmic time. With current in- 



tronomy, Inc. under the NASA contract NAS 5-26555. These observations 
are associated with programs 9051 and 12214. 
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strumentation, the most effective way to study this gas (in its neu- 
tral, ionised and molecular phases) is with absorption spectra to- 
wards bright, high redshift sources such as quasars and gamma- 
ray bursts (e.g. Rauch 1998; Vreeswijk et al. 2004; Wolfe, Ga- 
wiser & Prochaska 2005; Prochaska et al. 2007b; Ledoux et al. 
2009). The damped Lyman-a systems (DLAs) and sub-DLAs (with 
log Af(Hi) > 20.3 and 19.0 < log N(H i) < 20.3, respectively) have 
garnered particular attention as the galactic scale absorbers that 
contain the bulk of the neutral gas, and which can be readily iden- 
tified in relatively low resolution spectra (e.g. Wolfe et al. 1986; 
Storrie-Lombardi & Wolfe 2000; Ellison et al. 2001; Peroux et al. 
2003; Rao, Turnshek & Nestor 2006). There are now some 1000 
DLAs identified, with the largest haul coming from the Sloan Digi- 
tal Sky Survey (SDSS; Noterdaeme et al. 2009; Prochaska & Wolfe 
2009). High resolution spectroscopic follow-up has yielded chem- 
ical abundances for over 200 DLAs and sub-DLAs (e.g. Prochaska 
et al. 2007a; Akerman et al. 2005; Dessauges-Zavadsky, Ellison 
& Murphy 2009; Meiring et al. 2009). These abundance measure- 
ments have been brought to bear on several issues of nucleosyn- 
thetic origin and enrichment mechanisms (e.g. Pettini et al. 2000, 
2002; Ellison, Ryan & Prochaska 2001; Prochaska et al. 2002; 
Dessauges-Zavadsky et al. 2002; Ledoux, Bergeron & Petitjean 
2002; Cooke et al. 2011). 

Despite this wealth of spectroscopic data, surprisingly little is 
known about the optical emission properties (such as morphologies 
and luminosities) of DLA and sub-DLA galaxies. Whilst there has 
been some success in imaging the host galaxies at low redshifts 
(e.g. Le Brun et al. 1997; Chen & Lanzetta 2003; Chun et al. 2010; 
Rao et al. 2003, 2011), only a handful of high z galaxy counter- 
parts have been identified (e.g. Warren et al. 2001; Weatherley et 
al. 2005; Fynbo et al. 2010, 2011; Noterdaeme et al. 2012). To ad- 
dress the difficulty of detecting high redshift DLA galaxies, several 
novel observing techniques have been initiated, such as imaging be- 
low the Lyman limit of a higher redshift absorber or using integral 
field units (O'Meara, Chen & Kaplan 2006; Fumagalli et al. 2010; 
Peroux et al. 2011). At least at lower redshifts, where the sample 
sizes are sufficient to draw general conclusions, galaxies selected 
based on their H i properties alone seem to represent a wide vari- 
ety of galaxy morphologies and luminosities (e.g. Rao et al. 2003, 
2011). 

Ultimately, however, it is perhaps the physical conditions 
within the ISM that may reveal more about the nature and processes 
in high redshift galaxies. There are a few promising avenues for 
making such measurements in DLAs. For example, in cases where 
molecular hydrogen is detected (e.g. Ge & Bechtold 1997; Ledoux 
et al. 2003; Noterdaeme et al. 2008) the relative populations of the 
rotationally excited states can be used to constrain temperatures and 
densities in the ISM (e.g. Srianand et al. 2005; Noterdaeme et al. 
2007, 2010). In cases where the background QSO is radio-loud, an 
alternative measure of the temperature can be obtained through ob- 
servations of the Hi 21 cm line in absorption. In combination with 
the H i column density (measured from the Lya absorption profile), 
the Hi 21cm optical depth yields the so-called spin temperature 
(T,) of the neutral gas. For the general case of multiple phases of 
neutral gas along the line of sight (e.g. a cold phase with a kinetic 
temperature of « 40 - 200 K and a warm phase with a kinetic tem- 
perature of a 5000 - 8000 K; Wolfire et al. 2003), the derived spin 
temperature gives the column density weighted harmonic mean of 
the spin temperatures of the multiple phases, and thus the distribu- 
tion of neutral gas in different phases of the ISM (e.g. Kanekar & 
Briggs 2004). Since the detection of molecular hydrogen in DLAs 
is biased towards the highest metallicity systems (Petitjean et al. 



2006; Noterdaeme et al. 2008), the use of spin temperatures is an 
appealing option for a more general overview of conditions in the 
neutral ISM. 

Despite its promise, the collection of spin temperature mea- 
surements in DLAs has been fairly slow, with the first 25 years of 
Hi 21 cm absorption studies yielding only 10 detections in DLAs 
at all redshifts, with only 3 detections at z > 1.7 (e.g. Brown 
& Roberts 1973; Wolfe & Davis 1979; Wolfe, Briggs & Jauncey 
1981; Wolfe et al. 1985; Lane et al. 1998; Chengalur & Kanekar 
1999). The primary reason for the lack of spin temperatures is the 
paucity of DLAs at suitable redshifts in front of radio-loud QSOs. 
Moreover, the relatively narrow bandwidth of spectrometers has 
generally made blind H 1 21 cm absorption surveys impractical. The 
last decade has seen notable improvement in Hi 21 cm absorption 
statistics, due to progress on several fronts. First, optical spectro- 
scopic surveys that include large numbers of QSOs are now avail- 
able to cross-correlate with all-sky radio surveys. This approach 
has led to the detection of Hi 21cm absorption associated with 
DLAs and Mg n absorbers towards SDSS QSOs (e.g. Kanekar et 
al. 2009b; Gupta et al. 2009; Srianand et al. 2010, 2012). A more 
targeted approach is to design an optical search for DLAs towards 
radio-loud QSOs (Ellison et al. 2001; Jorgenson et al. 2006; Elli- 
son et al. 2008). To fully exploit these new surveys, there has also 
been progress on the technical side, with improvements in the fre- 
quency range accessible by large radio telescopes such as the West- 
erbork Synthesis Radio Telescope (WSRT), the Green Bank Tele- 
scope (GBT) and the Giant Metrewave Radio Telescope (GMRT). 
As a result of these advances, the opportunity now exists to gather 
significant numbers of T, measurements in DLAs at high redshifts, 
z>2. 

Almost immediately, the measurements of T 5 in high redshift 
DLAs yielded a surprise result: whereas z abs < 0.5 DLAs exhibit a 
range of spin temperatures, the sample of absorbers at z abs > 2 ap- 
peared to be dominated by spin temperatures in excess of 1000 K 
(Carilli et al. 1996; Chengalur & Kanekar 2000; Kanekar & Chen- 
galur 2003; Kanekar et al. 2007; Kanekar et al. 2012). Although 
York et al. (2007) discovered the first low T, DLA at z abs > 2, 
whether or not there is a wholesale evolution in T, with redshift has 
remained a hotly debated topic (e.g. Kanekar & Chengalur 2003; 
Wolfe, Gawiser & Prochaska 2003; Curran et al. 2005; Kanekar 
et al. 2009c). In addition to being limited by statistics, there has 
been considerable focus on the possible effect of the so-called cov- 
ering factor, / (e.g. Wolfe, Gawiser & Prochaska 2003; Curran et 
al. 2005; Kanekar et al. 2009a). If the background radio source is 
extended, then the foreground absorber will not fully cover it; this 
would imply a low covering factor, / << 1. Since most Hi 21 cm 
absorption studies actually determine the ratio TJf (and not T s 
alone), assuming / = 1 would result in an over-estimate of the 
spin temperature. It has also been claimed (Curran & Webb 2006) 
that due to cosmological geometry, different covering factors are to 
be expected in different redshift regimes. In particular, absorbers 
at high redshifts, z abs > 1.6, are always at larger angular diameter 
distances than the background QSOs, potentially leading to system- 
atically lower covering factors for higher redshift DLAs. However, 
in a Very Long Baseline Array (VLBA) study of radio-loud QSOs 
with foreground DLAs, Kanekar et al. (2009a) found none of the 
24 DLAs of their sample have / < 0.4, and that there was no dif- 
ference in the distributions of the covering factors of the high-z and 
low-z DLA samples. Kanekar et al. (2009a) hence concluded that 
low covering factors in high-z DLAs were not the cause of the ap- 
parent redshift evolution in the spin temperature. 

A second intriguing suggestion involving spin temperatures is 
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the possible anti-correlation between T, and metallicity (Kanekar 
et al. 2009c). Since metals are the primary coolant in the ISM, one 
would expect the fraction of cold gas to increase as the metallic- 
ity increases (Kanekar & Chengalur 2001). Kanekar et al. (2009c) 
reported a tentative detection of this anti-correlation at 3.6<r signifi- 
cance from a sample of 26 DLAs, of which 10 had measurements of 
(rather than limits on) both T, and metallicity. Again, the absorber 
covering factors were available for only a subset of the DLAs. 

The main challenge in confirming the redshift evolution of T s 
and the anti-correlation of T, and metallicity is simple statistics. For 
example, there remains a glaring gap in the redshift coverage of T s 
measurements between 0.6 < z a b s < 1-7, the so-called Hi 21 cm 
redshift desert. There are only two published Hi 21 cm detections 
in DLAs in this regime, which spans 4 gigayears in cosmic time 
(Brown & Roberts 1973; Kanekar et al. 2009b). To fill this regime, 
several surveys have searched for Hi 21cm absorption in strong 
Mgn systems that have been identified from ground-based opti- 
cal surveys (e.g. Ellison et al. 2004; Prochter, Prochaska & Buries 
2006). Although these Hi 21 cm surveys have now yielded a sig- 
nificant number of Hi 21 cm absorption detections in the redshift 
desert (Gupta et al. 2007, 2009; Srianand et al. 2008; Kanekar et al. 
2009b), the critical A^Hi) measurements required to determine the 
spin temperatures have been missing. A second obstacle has been 
determining covering factors for the full redshift range of H i 21 cm 
absorbers, since these corrections are required to confidently test 
any trends in T,. 

In this paper, we present a comprehensive set of multi- 
wavelength data for six Hi 21 cm absorbers in the redshift desert. 
This is an effective way to improve spin temperature statistics since 
the main limitation of identifying systems with Hi 21cm detec- 
tions has already been addressed. The data span six decades in rest 
wavelength ranging from near ultra-violet (UV) absorption at 1216 
A through the optical and out to the redshifted 21 cm line frequency 
in the radio regime. The sample is selected primarily based on 
Hi 21 cm absorption in a sample of strong Mgn absorbers (Gupta 
et al. 2009; Kanekar et al. 2009b). The near-UV data cover the Lya 
transition, from which the H i column density of the absorber can 
be measured. We also present VLBA images of all of the QSOs 
in order to determine their low-frequency covering factors . Finally, 
we complement the dataset with optical echelle spectra that permit 
the measurement of the column densities of a variety of elements. 
Combined with the N(Hi) measurements from the Lya data, we 
can hence also determine metal abundances. The dataset presented 
here is therefore unique at this redshift and represents the first sys- 
tematic study of metallicities and spin temperatures in the redshift 
desert. In this paper, we focus on the metallicities of H i 21 cm ab- 
sorbers. In a companion paper (Kanekar et al. in preparation) we 
will address the issue of T, evolution. 



2 OBSERVATIONS 



2.1 H i 21 cm spectroscopy with the GBT 

The targets in our sample are listed in Table[T] Of the 6 absorbers in- 
cluded in the current sample, 5 have previously published H 1 21 cm 
detections with the GMRT or the GBT (Kanekar et al. 2009b; Gupta 
et al. 2009). The final target, J1431+3952, is a new Hi 21 cm ab- 
sorption detection that we present here for the first time. For com- 



QSO 
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Zabs 


NUV 


Optical 


S 1.4 GHz 








mag. 


mag. 


(Jy) 


B0105-008 


1.374 


1.37078 


18.8 


g = 17.7 


0.93 


B0237-233 


2.223 


1.67235 


19.5 


B = 16.8 


6.26 


B0801+303 


1.451 


1.19110 


19.8 


g = 18.3 


1.27 


J1431+3952 


1.215 


0.60190 




8 = 16.5 


0.22 


J1623+0718 


1.648 


1.33567 




g = 17.7 


0.081 


B2355-106 


1.639 


1.17303 


20.7 


8 = 18.8 


0.77 



Table 1. The sample of QSOs with H i 2 1 cm absorption studied in this 
work. NUV magnitudes are taken from the GALEX survey, while the 
1.4 GHz flux densities are from the National Radio Astronomy Observa- 
tory (NRAO) Very Large Array (VLA) Sky Survey (NVSS; Condon et al. 
1998). All DLA absorption redshifts quoted in this paper are values derived 
from the peak H 1 21 cm absorption (Kanekar et al. 2009b; Gupta et al. 2009; 
see Table |2). 
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Figure 1. GBT spectrum of the H 1 2 1 cm absorption at Zubs = 0.602 towards 
J1431+3952. 



pleteness, the details of the Hi 21 cm detections of all six absorbers 
in our sample, including J1431+3952, are given in Table|2] 

We initially observed J143 1+3952 with the GBT on Septem- 
ber 29 2009, using the PF1-800 receiver, the Spectral Processor and 
a 2.5 MHz bandwidth centred at 886.645 MHz and sub-divided 
into 1024 channels. This was part of a large new GBT Hi 21 cm 
absorption survey of strong Mgn absorption systems (proposal 
AGBT10B_008; PI Kanekar). Position-switching with an On/Off 
time of 5 minutes was used to calibrate the bandpass, while the sys- 
tem temperature was determined via a blinking noise diode. The 
total on-source time was 15 minutes, sub-divided into 10-second 
integrations. 

The initial observations of J1431+3952 revealed a strong ab- 
sorption feature at the expected redshifted H 1 21 cm line frequency. 
We hence repeated the observations on October 20 2010 and Jan- 
uary 21 2011, using a 1.25 MHz bandwidth to obtain higher ve- 
locity resolution, as well as longer integration times (30 minutes 
on-source at each epoch). These observations confirmed the 21 cm 
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Figure 2. VLBA images of the 6 quasars in the sample, ordered (from top left) by quasar right ascension. The map frequency is at the top right of each panel. 



absorption by a detection of the doppler shift in the observed 
Hi 21 cm line frequency, due to the Earth's motion around the Sun. 

The GBT data were analysed using the NRAO AIPS++ 
single-dish package DISH, following standard procedures. The 
data were first visually inspected for radio frequency interference 
(RFI) and Spectral Processor failures, with corrupted data edited 
out. After editing, the individual spectra were calibrated and aver- 
aged together, in the barycentric frame, to measure the flux density 
of J1431+3952 at 886.64 MHz. A second-order spectral baseline 
was then fit to each 10-second spectrum, during calibration, and 
subtracted out, before averaging the residual spectra together, again 
in the barycentric frame, to produce the final spectrum. 

We note that very different source flux densities were obtained 
in the October and January observing runs (179 - 240 mjy), al- 
though the depth of the absorption features agree within the noise. 
The strong out-of-band RFI present above 887 MHz makes it diffi- 
cult to accurately calibrate the flux density scale in such total power 
measurements. We have used the average flux density from the two 
runs (210 mjy) to compute the Hi 21 cm optical depth; the final 
quoted error includes the above systematic error in the flux density 
scale calibration, assumed to be ±30 mjy, i.e. « 15%. 

The final H 1 21 cm spectrum from the higher-resolution obser- 
vations of 2010 October and 2011 January is shown in Fig. Q] this 
has a velocity resolution of 0.83 km s~' (after Hanning-smoothing 
and re-sampling) and a root mean square (RMS) noise of 4.5 mjy 



per independent 0.83 km s spectral channel. Two clear features 
are visible in the spectrum, with depths of ~ 40 mjy; each has 
high statistical significance and was detected on all three observing 
runs. The integrated H 1 21 cm optical depth is (3.07 ±0.34) km s _1 , 
where we have added in quadrature the errors on the integrated flux 
density and the flux density scale. 



2.2 High spatial resolution radio imaging with the VLBA 

We used the VLBA in proposals BK159 and BK174 (PI Kanekar) 
to obtain high spatial resolution images of five of the back- 
ground quasars of our sample, at frequencies close to the red- 
shifted Hi 21cm line frequency. The observations of four tar- 
gets (B0105-008, B0237-233, B0801+303 andB2355-106) used 
the 606 MHz receivers, while J1431+3952 was observed with the 
1.4 GHz receivers. The 606 MHz observations were carried out 
on 11 September 2009 and 4 January 2010 using two adjacent 
2 MHz baseband channel pairs, while the 1.4 GHz observations 
were carried out on 27 December 2011 using four adjacent 16 
MHz baseband channel pairs, both with right- and left-hand cir- 
cular polarizations, and sampled at 2 bits. The 606 MHz and the 
1.4 GHz data were correlated using the VLBA hardware correla- 
tor and the VLBA DiFX software correlator (Deller et al. 2011) 
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QSO 


Telescope 


Resolution 


On- source Time 


RMS 


JrdV 


Reference 






(km s -1 ) 


(hours) 


(mJy) 


(kms- 1 ) 




BO1O5-0O8 


GMRT 


2.0 


3.5 


3.3 


0.995 ± 0.023 


Kanekar et al. (2009b) 


B0237-233 


GBT 


0.69 


1.5 


56 


0.076 ± 0.016 


Kanekar et al. (2009b) 


B0801+303 


GMRT 


3.6 


5.0 


2.2 


0.305 ± 0.025 


Kanekar et al. (2009b) 


J143 1+3952 


GBT 


0.8 


1.0 


8.3 


3.07 ± 0.34 


This work 


J1623+0718 


GMRT 


3.9 


11.7 


1.0 


0.916 ± 0.099 


Gupta et al. (2009) 


B2355-106 


GMRT 


1.8 


4.0 


2.5 


0.256 ± 0.062 


Kanekar et al. (2009b) 



Table 2. Summary of H i 2 1 cm observations. 



QSO 


On-source 


Frequency 


Beam 


RMS 


S core 


S'wt 


Ang. size 


Spatial extent 


/ 




time (hours) 


(MHz) 


(mas 2 ) 


(mJy/Bm) 


(Jy) 


(Jy) 


mas X mas 


pc 2 




B0105-008 


2.0 


606 


45.3 x 18.1 


4.0 


0.403 


1.26 


12.2x0.0 


104x0 


0.32 


B0237-233 


2.0 


606 


33.8 x23.6 


23 


4.95 


5.53 


11.3 X9.7 


97 x 83 


0.90 


B0801+303 


2.0 


606 


29.6 x 12.0 


2.3 


0.053 


2.45 


12.1 X 4.7 


101 x39 


0.02 


J1431+3952 


1.6 


1420 


12.9 x 5.3 


0.14 


0.065 


0.207 


3.2x0.0 


21 xO 


0.32 


J1623+0718 


2.0 


1420 


18.4 X 6.3 


0.21 


0.028 


0.081 


4.3 x 1.4 


36 x 12 


0.34 


B2355-106 


2.0 


606 


44.4 X 22.8 


2.2 


0.42 


0.42 


10.8x2.2 


90 x 18 


1.0 



Table 3. Summary of VLBA observations and results. The covering factor, /, is the ratio of the flux in the core (Score) to me tota ' integrated flux (Sim). A 
spatial extent of zero indicates a size much smaller than the synthesised beam. 



in Socorro, NM, respectively, with a 4 second correlator integra- 
tion time and 32 spectral channels per respective baseband chan- 
nel. The total bandwidth of the data was 4 and 64 MHz for the 
606 MHz and 1.4 GHz observations, respectively. Phase referenc- 
ing was not used. Bandpass calibration was carried out using obser- 
vations of the strong calibrators 3C345, 3C454.3, 3C48 and 3C147. 
The total on-source time was ~ 2 hours for the 606 MHz targets and 
1.6 hours for J1431+3952. In addition, we also retrieved and anal- 
ysed an archival VLBA 1.4 GHz dataset on the last source in our 
sample, 11623+0718 (proposal BG0189, PI Gupta); these observa- 
tions had a similar setup to the above noted 1.4 GHz observations 
but used four adjacent 8 MHz baseband channel pairs, resulting in 
a total bandwidth of 32 MHz. 

The VLBA data were analysed in "classic" AIPS, using stan- 
dard procedures. After applying the initial ionospheric corrections 
and editing data affected by RFI, the flux density scale was cali- 
brated using the measured antenna gains and system temperatures. 
This was followed by fringe-fitting for the delay rates and band- 
pass calibration, after which the data towards each target were av- 
eraged in frequency to produce a single-channel visibility dataset. 
The 606 MHz data were found to be of relatively low quality, 
with only 6-7 working antennas and poor ionospheric conditions. A 
number of cycles of self-calibration and imaging were used to de- 
rive the antenna-based gains for each target, with phase-only self- 
calibration used for the 606 MHz datasets and amplitude-and-phase 
self-calibration for the 1 .4 GHz datasets. For each source, the above 
cycles were repeated until no improvement was seen on further 
self-calibration. The final images obtained from the above proce- 
dure are shown in Fig.ff] in order of increasing right ascension. For 
each source, the synthesised beam and off-source RMS noise are 
listed in Table[3] 

Finally, the AIPS task IMFIT was used to fit a gaussian model 
to the radio continuum images, to measure the compact flux den- 
sity of each quasar. A single gaussian was used as the model for 
all sources but J1431+3952 and B0105-008, where a 2-gaussian 
model was used due to the clear presence of a jet or extended struc- 
ture. The VLBA results are summarised in TablefS] whose columns 



contain (1) the source name, (2) the on-source time, (3) the observ- 
ing frequency (MHz), (4) the synthesised beam, (5) the RMS noise 
(mJy/Bm), (6) the integrated core flux density (Jy) from the VLBA 
images, obtained using JMFIT, (7) the total flux density from low 
resolution interferometric images at a similar frequency ( Jy), (8) the 
deconvolved angular size of the core radio emission, in mas x mas, 
(9) the spatial extent of the core emission at the absorber redshift, 
and (10) the covering factor /, obtained by taking the ratio of the 
core flux density to the total flux density. We do not quote formal 
errors on the covering factor for two reasons. First, the VLBA and 
low-frequency measurements were not carried out simultaneously, 
which means that it is hard to account for source variability in the 
error estimates. Second, it is possible that the foreground DLA also 
covers some of the extended radio emission, making the values of 
/lower limits. For two sources (B0105-008 and J143 1+3952), the 
deconvolved angular size and the spatial extent are 0.0 along one 
axis; this implies that the core is unresolved along this axis, with 
an angular size significantly smaller than the VLBA synthesised 
beam. 



2.3 Near TJV spectroscopy with HST/COS 

Hubble Space Telescope (HST) Cosmic Origins Spectrograph 
(COS) spectra have been obtained for 4 of the 6 absorbers in our 
sample (proposal ID 12214, PI Ellison). For B0237-233, the Lya 
line is covered by a Very Large Telescope (VLT) UV and Visual 
Echelle Spectrograph (UVES) spectrum and discussed in Section 
\2A\ The final target, J 143 1+3952, has an existing archival HST 
Space Telescope Imaging Spectrograph (STIS) spectrum from Cy- 
cle 10 (proposal ID 9051, PI Becker) for which we use the pipeline- 
reduced spectrum available from the archive. 

The HST/COS data were obtained with the near-UV (NUV) 
channel and the G230L grating which yields a full width at half 
maximum (FWHM) resolution of ~0.8 A. The target was acquired 
using Guide Star Catalog 2 (GSC2) coordinates in ACQ/IMAGE 
mode in order to minimise acquisition overheads. Spectroscopic 
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Q0105-008, z=1.3725, log N(HI) = 21.70±0.15 



J1431+3952 , z = 0.603, log N(HI) = 21.2±0.1 





1150 1200 1250 

Wavelength (A) 
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Figure 3. Voigt profile models to the hya line for the absorbers in our sample are shown as red curves. The dashed lines indicate the uncertainties given in the 
panel captions. 



data were obtained in TIME-TAG mode with at least two exposures 
per target, each in a different FP-POS position. Total exposure times 
and central wavelength settings are given in Table |4] 

The COS spectra were homogeneously reduced using CAL- 
COS v2.15.4 and custom routines for background estimation 
and co-addition of individual exposures in the Poisson limit 
(e.g. Worseck et al. 2011). The data were flat-fielded, deadtime- 
corrected, doppler-corrected and wavelength-calibrated using the 
default CALCOS calibration files as of September 2011. For 1-d 
extraction we adjusted the source extraction window to a rectangu- 
lar box (height 17 pixels, corresponding to 5x the spatial FWHM) 
to preserve integer extracted counts in the Poisson regime, and 
to maximise the signal-to-noise ratio (S/N) for our well-centred 
point sources. Likewise, the background extraction windows were 
maximised (150 pixels each), and background smoothing by the 
pipeline turned off to investigate spectral variations of the back- 
ground. After visual inspection, the smoothed background was de- 
termined by a 51 -pixel running average in wavelength. Individual 
exposures were co-added by adding the integer source counts and 
the non-integer mean background per pixel on the homogeneous 
NUV wavelength grid, accounting for individual pixel exposure 



times induced by offsets in dispersion direction. The counts were 
then flux-calibrated with the NUV flux calibration curve and the 
pixel exposure times. The S/N ratios per pixel (in the continuum on 
either side of the hya absorption) of the final combined spectra are 
given in Table [4] 

The Starlink package DIPSO was used to model a Voigt- 
profile absoiption to the normalised spectra. To normalise the spec- 
tra, a cubic spline was fit through unabsorbed regions of the QSO 
continuum and the data then divided by the polynomial fit. Despite 
the relatively low S/N, the A'(Hi) could usually be determined to 
within + 0.1 dex, which is typical for moderate resolution spec- 
tra of low S/N (e.g. Russell, Ellison & Benn 2006). The dominant 
source of error in N(H i) is from uncertainty in the continuum place- 
ment. This is a particular issue for the proximate DLA towards 
B0105-008 where the Lyo- absorption lies on top of the QSO's 
Lya emission. Nonetheless, errors on A^(H i) are less < 0. 15 dex for 
all the bona-fide DLAs (log iV(Hi) > 20.3) in our sample. Tableg] 
lists the measured A?(Hi) column densities for both the COS and 
STIS data and Figure[3]presents the Lya fits to the absorbers in our 
sample. 

From Table [4] it can be seen that 4/6 absorbers in our sample 
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qualify as DLAs; indeed these 4 have high Hi column densities, 
in excess of 10 21 citT 2 . The absorbers towards B0801+303 and 
B0237-233 have 7V(Hi) values below the DLA threshold, but can 
still be well fit with reasonable accuracy and fall into the category 
of sub-DLAs. 



QSO 


Zabs 


log N(HI) 


T s (K) 


[Zn/H] 


B0105-008 


1.37078 


21.70±0.15 


305±45 


-1.40 


B0237-233 


1.67235 


19.78±0.07 


390±127 


-0.57 


J1431+3952 


0.60190 


21.2+0.1 


90±23 


-0.80 


J1623+0718 


1.33567 


21.35+0.10 


460±105 


-1.07 


B2355-106 


1.17303 


21.0±0.1 


2145±570 


-0.87 



2.4 Optical echelle spectroscopy with Keck/HIRES and 
VLT/UVES 

High resolution echelle spectra have been obtained for B0105-008, 
B1430-178, J1431+3952 and J1623+0718 using the High Res- 
olution Echelle Spectrograph (HIRES; Vogt et al. 1994) on the 
Keck telescope and UVES on the VLT. Our HIRES observations 
of B0105-008 were supplemented with archival UVES data of this 
target (proposal 082.A-0569, PI Srianand) in order to extend the 
wavelength coverage of the spectrum. Archival UVES data were 
also downloaded and reduced for B2355-106 (proposal 085.A- 
0258, PI Srianand). The same UVES program also obtained data 
for J1623+0718 which we use to supplement our HIRES data 
(which extend further into the blue and cover the UVES inter-arm 
and chip gaps). B0237-233 was a target in the ESO Large Program 
166.A-0106 (PI Bergeron) and its UVES data cover both the Lya 
and metal transitions. We do not have echelle spectra for the final 
target in our sample (B0801+303). Table [5] summarises the avail- 
able observations, instrument set-ups, wavelength coverage and fi- 
nal S/N ratios. 



UVES is a dual-arm echelle spectrograph with a grating cross 
disperser (Dekker et al. 2000). The data presented in this paper 
were obtained in service mode (or from the archive) using various 
combinations of standard wavelength settings. The actual wave- 
length coverage is given in Table [5] as measured directly from the 
spectra. The CCDs were binned 2x2 and a 1 arcsecond slit was 
used, resulting in a typical resolution of R = 42,000. The spec- 
tra were reduced using a custom version of the Midas reduction 
pipeline, with manual intervention at several stages. In particular, 
the tolerance of the wavelength solution was adjusted to optimise 
the number of identified arc lines and accuracy of the wavelength 
solution. The typical accuracy of the wavelength solution was 60- 
70 m/s. The extracted spectra from each exposure were converted 
to a common vacuum-heliocentric wavelength scale. The heliocen- 
tric correction was calculated externally, as the UVES fits headers 
do not always accurately record this parameter. Where multiple ex- 
posures in a given setting were obtained, or where there is spectral 
overlap between settings, the data were combined, weighting by 
the inverse of the flux variances. The final extracted 1-dimensional 
spectrum was normalised by fitting the continuum with Starlink's 
DIPSO package. 

HIRES data for three targets were obtained in PI mode. 
B0 105 -008 was observed with HIRES on 12 January, 2008 with 
the blue cross-disperser and collimator (i.e. HIRESb). We em- 
ployed the CI decker for a nominal spectral resolution of FWHM ~ 
6 km s -1 . Two exposures totalling 7200s were acquired under fair 
conditions. J1431+3952 was observed on 16 January 2012 in good 
conditions, also with the CI decker and HIRESb, with a single ex- 
posure of 1200s. The HIRES data were reduced with the HIRedux 
pipeline that is available within the XIDL software packagfl The 

1 http://www.ucolick.org/~xavier/IDL 



Table 8. Summary of DLA properties. Spin temperature values take into ac- 
count the covering fractions determined from VLBA mapping (see Section 

£2). 



main difference in the processing of the UVES and HIRES data 
occurs in the normalization procedure. Whereas normalisation is 
the final step in the UVES reductions, the HIRedux package fits 
the continuum order-by-order before conversion to a 1-dimensional 
spectrum. The HIRES spectra are also corrected to a vacuum helio- 
centric scale. 

Mosaics of many of the metal line transitions detected in the 
HIRES and UVES data are presented in Figures lA"Tl to lA5l 



2.5 Metallicities and spin temperatures 

Metal column densities were determined from the echelle spectra 
using the apparent optical depth method on unsaturated transitions 
(Table[6]l. The metal column densities were combined with /V(Hi) 
following standard procedures to determine an abundance relative 
to solar, with solar abundances taken from Asplund et al. (2009): 



[xi m = \ og N(x) - log mm - pogm - io g (/ry . 



(i) 



We follow the standard assumption that each element's ab- 
sorption is dominated by a single ionization state such that N(Fe) 
= N(Fen) etc. This should be a good apprioximation for DLA gas 
shielded from the UV background and whose ionisation state is de- 
termined by its stellar population. Final abundances are given in 
Table UJ 

For H i 21 cm absorption studies of DLAs towards radio-loud 
QSOs, the Hi 21 cm optical depth, Hi column density /V(Hi) and 
spin temperature T, are related by the expression 



N(HI) = 1.823 x 10 1 



rdV, 



(2) 



where /V(Hi) is in cm , T s is in K, and dV is in km s . We 
adopt the measurements of A'(Hi), r and / derived in Sections [2.H - 
12.41 in the final estimates of our spin temperatures. In Table [8] we 
provide the final spin temperatures (corrected for the source cov- 
ering factor) for five of the targets of our sample. Note that the 
table does not list a spin temperature for the sub-DLA towards 
B0801+303; see Section 12.6.31 below. In Figure|4]we overlay the 
Hi 21 cm and metal line profiles. In all cases, we show the same 
two metal line transitions: Fe n /12260 and Mg i A 2852. Whilst Mg i 
is likely to be a better tracer of the cold neutral medium (CNM) gas 
that is associated with Hi 21 cm absorption than the singly ionised 
species that are the most commonly detected (e.g. Mg n, Fen etc.), 
it can also have contributions from warm ionised gas (e.g. Kanekar 
et al. 2010). However, in the absence of, for example, C i detections, 
Mg i is our best tracer of the CNM amongst the UV resonance lines. 
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QSO 


Telescope/ 


Grating/ 


# 


Time 


S/N 


log N(HI) 




Instrument 


wavelength (A) 


orbits 


(s) 


pix 




BO1O5-0O8 


HST/COS 


G230L/2950 


1 


2,454 


4 


21.70±0.15 


B0237-233 


VLT/UVES 


346+580 




42,927 


40 


19.78±0.07 


B080I+303 


HST/COS 


G230L/2635 


2 


5,619 


5 


19.0±0.3 


J1431+3952 


HST/STIS 


G230L 


1 


2,292 


15 


21.2+0.1 


J1623+0718 


HST/COS 


G230L/2950 


4 


11,499 


5 


21.35±0.10 


B2355-106 


HST/COS 


G230L/2635 


4 


11,451 


4 


21.0±0.1 



Table 4. Summary of Lya observations. The S/N is measured in the continuum on either side of the Lya absorption and is derived from the RMS noise. 
However, at low S/N (the Poission limit) the formal S/N is about a factor of two lower. 



QSO 


Instrument 


Date/ 


Time 


Instrument 


Wavelength 


S/N 






Archive ID 


(s) 


set-up 


coverage (nm) 


pix -1 


B0105-008 


HIRES 


Jan. 2008 


7,200 


0/1.250 


307-593 


10-20 




UVES 


082.A-0569 


14,790 


390+580 


350-450, 476-684 


20-50 


B0237-233 


UVES 


166.A-0106 


42,927 


346+580 


305-1008 


40-150 








48,030 


437+860 


(combined) 




J1431+3952 


HIRES 


Jan. 2012 


1,200 


0/1.156 


317-603 


6-15 


J1623+0718 


HIRES 


Sept. 2010 


7,200 


0/1.0005 


320-588 


10 




UVES 


085.A-0258 


13,356 


390+580 


349-451,479-680 


10-40 


B2355-106 


UVES 


085.A-0258 


23,373 


390+580 


350-450, 476-684 


8-25 



Table 5. Summary of our optical echelle observations. The quoted S/N range is over the full spectral coverage. For spectra obtained directly by us, we provide 
the date of the observations. For UVES archival data, the proposal ID is given. The instrument set-up for HIRES refers to the echelle/cross disperser angles. 



Below, we briefly summarise the noteworthy features of each 
absorber. 

2.6 Notes on individual absorbers 

The H 1 21 cm absorbers presented in this paper are part of two sam- 
ples of "strong" Mg n absorbers, selected from the SDSS and the 
CORALS survey (Ellison et al. 2001), towards radio-loud quasars 
that were surveyed for Hi 21 cm absoiption with the GMRT and 
the GBT (Kanekar et al. 2009b; Gupta et al. 2009). The selection 
criteria were (1) W^ 2196 > 0.5 A and W* 2600 > 0.34 A (where 
W* 2196 and IVq 12600 are the rest equivalent widths in the Mg ni2796 
and Fe H/12600 transitions, respectively) and (2) that the absorber 
should lie towards a radio-loud background quasar with flux den- 
sity > 300 mjy at the redshifted H i 21 cm line frequency. 

2.6.1 B0105-008 

The relative velocity between the DLA and the background QSO is 
Av ~ 200 km s _1 . Although this velocity offset is likely to be uncer- 
tain by several hundreds of km s -1 (van den Berk 2001; Richards et 
al. 2002; Shen et al. 2007), the absorber is almost certainly a "proxi- 
mate" DLA (PDLA; e.g. the surveys by Ellison et al. 2002; Russell, 
Ellison & Benn 2006; Prochaska, Hennawi & Herbert-Fort 2008). 
There is evidence that high N(H i) PDLAs may have elevated abun- 
dances, by a factor ~ 3 (Ellison et al. 2010, 201 1). Further, the spin 
temperatures in PDLAs are likely to be affected by the proximity 
of the gas to a bright radio source (e.g. Field 1958). We discuss the 
treatment of PDLAs in our sample explicitly in Section [3~2l 

The Hi 21 cm absorption towards this absorber has a narrow 
profile that is well fit with a simple 2-gaussian model (Gupta et al. 
2007; Kanekar et al. 2009b). As noted in Section [2~2l there is evi- 
dence of a core-jet structure in the 606 MHz VLBA image, which 
leads to a core covering factor of 0.32. The core is unresolved along 



the declination axis, with angular extent << 18 mas (i.e. spatial ex- 
tent « 150 pc). 

The Hi column density in this absorber has a larger uncer- 
tainty than that of the other DLAs in our sample as the absorption 
overlaps with the Lya emission line; this is a common problem for 
PDLAs. Nonetheless, the /V(Hi) is clearly extremely high in this 
absorber, the largest in our sample with log /V(Hi) = 21.70 + 0.15. 

Both HIRES and UVES data are available for this sightline; 
the latter were obtained from the UVES archive, from proposal 
082.A-0569 (PI Srianand). The UVES data have a higher S/N 
and wider wavelength coverage and are used in preference to the 
HIRES data for most transitions. However, the HIRES data have 
more extensive coverage of the blue spectrum for transitions with 
rest wavelengths 1290 < A <1475 A. 

The combination of good S/N and high N(H i) facilitates the 
detection of a wide range of metal line transitions, which all show 
a blended two-component structure. The Hi 21cm absorption is 
offset from the strongest metal line component, and lies between 
the two components in velocity space (Figure|4j. The column den- 
sity of Si ii is reported as a lower limit due to the mild saturation 
of the A 1808 transition. Overall, this absorber has a metallicity of 
Z ~ 1/25Z Q , which is relatively metal-poor for this redshift (see 
Section[3). The low metallicity also puts this absorber at odds with 
the majority of other high N(H i) PDLAs, which tend to have rela- 
tively high abundances (Ellison et al. 2010, 201 1). 

2.6.2 B0237-233 

The UVES spectrum of B0237-233, obtained as part of a UVES 
Large Program, (166.A-0106, PI Bergeron), is truly exquisite. 
The combination of very long integration times (approximately 25 
hours split between two different wavelength settings) and a bright 
apparent magnitude (B = 16.8) leads to S/N ratios up to 150. The 
high quality of the UVES spectrum, which extends as blue as 3050 



© 0000 RAS, MNRAS 000, 000-000 



DLAs in the redshift desert 



QSO 


-abs 


N(FeII) 


N(Sffl) 


N(ZnII) 


N(CrII) 


N(MnII) 


N(TiII) 


N(NiII) 


B0105-008 


1.37078 


15.59±0.03 


> 15.85 


12.93±0.04 


13.92±0.03 


13.29+0.04 


13.0+0.2 


14.25+0.04 


B0237-233 


1.67325 


14.57±0.02 


14.92±0.02 


11.84±0.09 


12.72±0.06 


12.24+0.06 


<11.46 


13.51+0.04 


J143 1+3952 


0.60190 


15.15±0.11 




13.03+0.19 


13.53+0.18 


13.07+0.14 


12.70+0.15 




J1623+0718 


1.33567 


15.28±0.05 


15.78±0.05 


12.91±0.09 


13.55+0.06 


13.06+0.06 


<12.54 


14.11+0.07 


B2355-106 


1.17303 


15.08±0.10 


15.42±0.13 


12.76±0.17 


<13.2 


12.91+0.10 


<12.9 


14.17+0.12 



Table 6. Column densities derived from the apparent optical depth method. No high resolution spectrum is available for B0801+303. For non-detections, 
quote 3<r upper limits. Lower limits are derived from lines exhibiting some degree of saturation. 



QSO 


£abs 


[Fe/H] 


[Si/H] 


[Zn/H] 


[Cr/H] 


[Mn/H] 


[Ti/H] 


[Ni/H] 


B0105-008 


1.37078 


-1.56 


> -1.36 


-1.40 


-1.42 


-1.89 


-1.61 


-1.65 


B0237-233 


1.67325 


-0.66 


-0.37 


-0.57 


-0.70 


-1.02 


< -1.23 


-0.47 


J1431+3952 


0.60190 


-1.50 




-0.80 


-1.31 


-1.61 


-1.41 




J1623+0718 


1.33567 


-1.52 


-1.08 


-1.07 


-1.44 


-1.77 


< -1.72 


-1.44 


B2355-106 


1.17303 


-1.37 


-1.09 


-0.87 


< -1.44 


-1.57 


< -1.01 


-1.03 



Table 7. Elemental abundances. No high resolution spectrum is available for B0801+303. 
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Figure 4. Velocity profiles for the five absorbers in our sample for which we have both Hi 21 cm and echelle spectra. In each panel, we show the sa 
transitions on the velocity scales used in Figures IAT1 - IA5 1 The Hi 21 cm profiles are arbitrarily offset and scaled for display purposes. 
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A, therefore permits a precise measure of A^H i) in this low z ab- 
sorber without the need for HST. We note the presence of a satu- 
rated Lya line in the red wing of the DLA. Lack of coverage of 
unsaturated higher order Lyman lines precludes an accurate A'(Hi) 
measurement of this line. However, this contamination is unlikely 
to affect the DLA fit, which is driven by the extended damping 
wings. 

Despite the low A'(Hi) of this absorber, which is below the 
traditional threshold of DLAs, the metal lines are sufficiently strong 
to enable us to determine a wide suite of abundances. The only 
species listed in Table [7] that we do not significantly detect is the 
challenging Ti n. There is actually a very weak absorption feature 
at the expected wavelength of Tin A 1910.6, but its significance is 
only 2a. We therefore conservatively report the formal 3cr limit for 
this transition. 

The Hi 21cm absorption is well aligned with the strongest 
of the two metal line components (Figure l4pl . There is a single 
pixel feature in the H i 21 cm absorption spectrum which is aligned 
with the weaker metal feature. Although suggestive, this feature is 
not statistically significant in the current GBT spectrum and is not 
included in the integrated H i 21 cm optical depth. 

The VLBA map of this quasar reveals a compact radio core, 
from which we derive a very high covering factor of 0.90. The spa- 
tial extent of the core emission is ~ 97 x 83 pc 2 , at the absorber 
redshift. 

2.6.3 B0801+303 

Kanekar et al. (2009b) discuss the likely low covering fraction of 
this absorber (see also Kunert et al. 2002). Our VLBA imaging 
shows a weak, compact continuum source with a flux density that 
is only about 2% of the total flux density measured with a lower 
resolution interferometer. This clearly suggests that the radio emis- 
sion from this source is extended on scales larger than 0.22 arc- 
sec; a value set by the VLBA short spacing limit at 606 MHz. It 
is therefore likely that the Hi 21 cm absorption arises against the 
extended radio flux and not against the radio core. The likely non- 
coincidence of the optical and radio sightlines implies that we can- 
not determine the spin temperature of this absorber. We therefore 
do not include this sub-DLA in our spin temperature analysis (it is 
excluded from Table [8} and have not sought to obtain high resolu- 
tion optical spectroscopy. 

2.6.4 J1431+3952 

This paper presents a new detection of Hi 21cm absorption in 
the sightline towards J1431+3952. The integrated Hi 21 cm opti- 
cal depth is quite high, with the profile consisting of 2 narrow ab- 
sorption components separated by ~ 15 km s -1 (Figure[T). The two 
Hi 21 cm absorption components are well aligned with the veloc- 
ity structure seen in Mg i (Figure|4j. The relatively low S/N around 
the Fe n A 2260 line makes its velocity structure harder to discern. 
However, the two component structure is still present in the singly 
ionised species, e.g. in Znn A 2026, again well aligned with the 
Mgi. 

The VLBA image of J143 1+3952 shows a core-jet structure, 
with a clear southern jet. The central radio emission was hence fit 
with a 2-component gaussian model and a covering factor of 0.32 

2 The extended absorption around the Mg i profile is due to contamination 
from the atmospheric A band. 



was obtained for the unresolved core. The core is extremely com- 
pact, with a deconvolved angular size of mas along the declination 
axis, indicating that the angular extent along this axis is << 5.3 mas 
(i.e. a spatial extent of « 35 pc). 

After correction for its sub-unity (~ 30%) covering factor, this 
DLA has the lowest measured spin temperature to date, T s = 90 ± 
23 K. The only other DLAs with spin temperatures < 150 K (after 
correcting for the covering factor) are the systems at z = 0.395 to- 
wards Q1229-021 (Brown & Spencer 1979) and z = 2.289 towards 
TXS031 1+430 (York et al. 2007) which have T,=95 and 120 K re- 
spectively. In the measurements of Galactic spin temperatures by 
Kanekar et al. (2011), two Milky Way sightlines have comparably 
low values (T,=89 and 1 19 K), both at very low Galactic latitudes. 
Given the suggestion that DLAs with low spin temperatures are of- 
ten (at least at low z) associated with spiral galaxies (Chengalur & 
Kanekar 2000; Kanekar & Chengalur 2001), the high implied CNM 
fraction towards J143 1+3952 may be associated with the intersec- 
tion of a galactic disk. 

The SDSS image of J1431+3952 shows a faint galaxy at ~ 
5 arcsec separation from the QSO. At a redshift z = 0.602 this 
would correspond to a physical projection of ~ 38 kpc, which is 
fairly large given the implied intersection of a cold disk gas. More- 
over, a fit of the SDSS photometry yields a photometric redshift 
z = 0.08 ± 0.02 (Michael Palmer, 2012, private communication). A 
high spatial resolution image of the J143 1+3952 field would be of 
great interest to attempt to identify the absorbing galaxy. 

2.6.5 J 1623+07 18 

The redshift of the Hi 21 cm absorption (z = 1.33567) is in excel- 
lent agreement (~ 4 km s ) with the centroid of the metal lines 
(z = 1.33570). For the majority of metal line transitions, we use the 
UVES spectrum, for its superior S/N. However, the Zn n A 2026 line 
falls between the blue and red arm coverage, so the HIRES data is 
used for this transition. The overall metallicity of this absorber is 
Z ~ l/10Z o . 

In order to compare the structure in the metal lines with the 
Hi 21 cm absorpion (originally presented by Gupta et al. 2009), we 
have re-reduced the archival GMRT data (project 14RSA01, PI Sri- 
anand) for this target using standard procedures in AIPS (see, for 
example, Kanekar et al. 2009b for a full description of the anal- 
ysis procedure of 610-MHz GMRT data). Our re-reduction of the 
GMRT data is presented in Figure [4] the strongest Hi 21 cm ab- 
sorption is well-aligned with the strongest metal component. 

The VLBA image of J 1623 +07 18 shows a strong core, with 
a possible weak jet, only tentatively detected in the present map. 
Only one-third of the total radio emission stems from the core, with 
a covering factor of ~ 0.34; most of the emission must arise from 
extended structure that is not detected in the VLBA image. The 
spatial extent of the core emission is quite small, ~ 36 x 12 pc 2 , at 
the absorber redshift. 



2.6.6 B2355-106 

The DLA towards B2355-106 has T,=2145 K, making it one of 
the highest spin temperatures to date in a DLA with a detection 
of Hi 21cm absorption. The Hi 21cm absorption reported by 
Kanekar et al. (2009b) in this DLA is located in a single narrow 
component at z = 1.17303. In contrast, the metals show absorption 
over ~ 150 km s~' (Figure lA5l . with the strongest Fen absorption 
at 1.17230, offset from the Hi 21 cm absorption by ~ 100 km s -1 , 
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Figure 5. Spin temperature as a function of redshift for all sub-DLAs 
and DLAs (including proximate absorbers) with measured covering factors. 
Filled and open circles represent H i 21 cm detections and limits from the 
literature respectively. Our intermediate redshift sample is shown by filled 
stars in the redshift desert (0.6 < z < 1.7, shaded region). Red stars indicate 
DLAs and blue stars show sub-DLAs. The horizontal dashed line shows a 
spin temperature of 400 K. 



However, Mgi is strongest in the redder components (at v > 
km s" 1 in Figure IaBI . which matches well the Hi 21 cm absorp- 
tion (Figure[4]l. Although not statistically significant, there may be 
some extended H 1 21 cm absorption bluewards of the main detected 
feature, spanning the broad velocity range over which the metal ab- 
sorption is observed (Figure |4). However, the extended absorption 
is not detected in the independent spectrum of Gupta et al. (2009), 
so we do not include it in the integrated optical depth and spin tem- 
perature determination. 

There is possible Cr nA 2056 absorption, but the significance 
of the feature is less than 3 cr, so a limit is quoted for N(Crn) in 
Table|7] The Tin /11910 line is formally detected at 4.9 cr, but vi- 
sually this is not very convincing. The significance of detections 
in this absorber are complicated by the broad velocity structure, so 
we conservatively quote the Ti n as a limit. The same broad struc- 
ture, combined with a relatively low S/N, makes weak features par- 
ticularly uncertain, which may lead to over-estimation of column 
densities (e.g. Zn n). 

The VLBA image of this source shows an extremely compact 
radio core, with a spatial extent of ~ 90 x 18 pc 2 at the absorber 
redshift, and a covering factor of unity. 



3 RESULTS AND DISCUSSION 
3.1 Filling in the redshift desert 

The measurements presented in this paper are combined with a 
wider literature sample of T , measurements in order to provide the 
first attempt at bridging the high and low redshift observations and 
to improve the overall statistics of spin temperatures and metallic- 
ities. The spin temperatures of DLAs and sub-DLAs are shown as 
a function of redshift in Figure [5] All spin temperatures in this fig- 
ure include the estimates of the DLA covering factor; indeed, in the 
following sections, we will deal only with /-corrected spin temper- 



atures unless otherwise statecJJ In a companion paper (Kanekar et 
al. in preparation), we will present the full DLA sample, including 
both new GBT/GMRT Hi 21cm absorption studies and sources 
from the literature, as well as the results of an extensive VLBA 
campaign to determine covering factors for many of the absorbers 
presented in Figure[5] The redshift desert between 0.6 < Za>s < 1.7 
is obvious in Figure [5]- prior to the work presented here, only 2 
Hi 21 cm detections existed in this range. We have added a fur- 
ther 5 T s measurements (filled stars) which sample the regime over 
which the spin temperature has been claimed to evolve. It is also 
interesting to note that our sample encompasses absorbers with a 
wide range of spin temperatures, including the lowest, and one of 
the highest values amongst Hi 21 cm detections. The full analysis 
of the spin temperature evolution will be presented by Kanekar et 
al. (in preparation). For the remainder of this paper, we focus on 
the metallicity properties of Hi 21 cm absorbers. 

3.2 Proximate DLAs and sub-DLAs 

Four of the DLAs shown in Figure [5] are proximate DLAs within 
3000 km s- 1 of the QSO (B0105-008, Q0405-331, Q1013+615 
and Q0528-250), including one of the absorbers in our intermedi- 
ate redshift sample (B0105-008). There are various reasons why 
PDLAs might be excluded from a general analysis of abundances 
in Hi 21 cm aborbers. First, radiation from the QSO may alter the 
hyperfine level populations in the Hi atoms (e.g. Field 1958; Wolfe 
& Burbidge 1975). The result of such irradiation would be that the 
observed spin temperature would depend both on the radiation field 
and on the distribution of gas in different temperature phases, sig- 
nificantly complicating the interpretation. Three of the four PDLAs 
in our sample have T,> 750 K; they are all at z abs > 2.5 where 
such spin temperatures are typical (Figure |3J. The final PDLA is 
the absorber towards Q0105-008 from our redshift desert sample 
(Zabs = 1-37), with T, = 305 K (including the correction for a cov- 
ering factor of 0.32). 

Independent of the H 1 21 cm absorption, it has been suggested 
that PDLAs with log Af(H i) > 2 1 are biased towards higher metal- 
licities (Ellison et al. 2010, 2011). Combined with a possible anti- 
correlation between T , and metallicity, the high N(H i) PDLAs may 
have preferentially low values of T,. Two of the PDLAs in the sam- 
ple, B0105-008 and Q0528-250, have log iV(Hi) > 21. The abun- 
dances of B0105-008 are presented here, and we find a very typi- 
cal metallicity for this redshift, [Zn/H]=-1.40. Q0528-250 was in- 
cluded in the PDLA study of Ellison et al. (2010) and does show a 
relatively high metallicity for its N(H i) (see Figure 15 in Ellison et 
al. 2010). 

In view of these potential biases, we limit the following anal- 
yses to absorbers that are at least 3000 km s from the published 
quasar systemic redshift. 

Another class of absorbers that is often distinguished from 
the general DLA population consists of the so-called sub-DLAs, 
systems with 19.0 < log iV(Hi) < 20.3. Although their AT(Hi) 
falls below the classical DLA criterion, they show clear damping 
wings and in many cases have reliable abundances (e.g. Dessauges- 
Zavadsky et al. 2003, although see Richter et al. 2005; Quast et al. 
2008; Milutinovic et al. 2010 for counter-examples). Although Per- 
oux et al. (2003) have suggested that sub-DLAs make an important 
contribution to cosmic Hi, particularly at high redshifts, O'Meara 

3 We use only covering factors that have been determined at frequencies 
within a factor of ~ 2 of the H i 21 cm absorption 
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et al. (2007) calculated that systems with log A'(Hi) < 20.3 repre- 
sent only ~ 10% of cosmic neutral hydrogen. Prochaska, Herbert- 
Fort & Wolfe (2005) have also argued that the sub-DLAs are not 
tracers of the main gas reservoirs for cosmic star formation. In- 
deed, there is evidence from Galactic sightlines that below the 
DLA Af(Hi) threshold there is only a small contribution from the 
CNM (Kaneker et al. 201 1). Finally, it is observed that low-redshift 
sub-DLAs (which have been mostly identified through their large 
Mgn equivalent widths) have systematically higher metallicities 
than DLAs at the same redshift (e.g. Peroux et al. 2008; Meiring 
et al. 2009; Dessauges-Zavadsky et al. 2009). 

Despite these literature debates, we do include sub-DLAs in 
the current work, as it is of interest to see whether the Galactic 
trends with A^H i) are also evident in cosmological DLAs and sub- 
DLAs. In particular, if the CNM fractions (and hence, spin tem- 
peratures) depend on metallicity, we might expect trends between 
these parameters to extend into the sub-DLA regime, regardless of 
the nature of the intervening galaxy. However, we do distinguish 
sub-DLAs from DLAs in figures and quoted statistics. 



3.3 Metallicities in H i 21 cm absorbers 

In Figure [6] we show the abundances of three commonly measured 
elements (Zn, Fe and Si) in DLAs and sub-DLAs with Hi 21 cm ab- 
sorption studies, as a function of redshift. We distinguish low and 
high (including lower limits) spin temperatures by plotting filled 
stars for absorbers with T, < 500 K and open stars for absorbers 
with higher values of T,. For comparison, we show a literature com- 
pilation (Berg et al, in preparation; Rafelski et al. 2012) from DLAs 
and sub-DLAs (open circles) towards all quasars (i.e. not merely ra- 
dio quasars). DLAs and sub-DLAs for both samples are shown in 
red and blue, respectively. 

The now well-documented (e.g. Peroux et al. 2008; Meiring et 
al. 2009) enhancement in sub-DLA abundances is clearly present 
at low redshifts (z abs < 1.7). Whilst it has been suggested that their 
high metallicities indicate that sub-DLAs sample the most massive 
galaxies, Dessauges-Zavadsky et al. (2009) have argued that this 
may be a selection effect, since the majority of sub-DLAs in this 
redshift range were selected via the Mg n-based survey of Rao et 
al. (2006). 

Whilst Mg n selection seems to sample DLAs in an unbiased 
way, Dessauges-Zavadsky et al. (2009) have argued that at lower 
yV(Hi) there may be a bias towards larger velocity widths and/or 
metallicities (see also Rao et al. 2006 for a discussion on selec- 
tion bias, and Bouche 2008 for further evidence of a bias towards 
high metallicity sub-DLAs selected by Mgn). In either case, the 
sample selected using Hi 21cm absorption (many of which, at 
Zabs < 1-7, were also originally selected based on Mgn rest equiv- 
alent width) show evidence of the same metallicity trends as the 
absorbers towards the full quasar sample. The H i 21 cm absorber 
metallicities cover the same ranges as those of the full sample, and 
the Hi 21 cm sub-DLAs show evidence for slightly higher metal- 
licities than DLAs at the same redshift. 

Although the elevated metallicities in Hi 21 cm sub-DLAs 
are consistent with optical samples, known trends in Hi 21cm 
absorbers may make the high metallicities in sub-DLAs with 
Hi 21cm absorption seem puzzling. Kanekar et al. (2011) have 
shown that Galactic sightlines with N(H i) in the sub-DLA regime 
tend to have high spin temperatures, and Kanekar et al. (2009c) find 
that high spin temperatures are more usually found in low metal- 
licity systems. So, we might expect that Hi 21 cm absorbers in the 



sub-DLA regime should have relatively low (rather than the ob- 
served high) abundances. 

In Figure|7]we put these arguments on a quantitative footing 
by showing the spin temperature and A?(Hi) Galactic data points 
used by Kanekar et al. (2011), as well as the values for Hi 21 cm 
detections in DLAs. The latter have been colour-coded by metal- 
licity. Metallicities are taken to be, in order of preference, [Zn/H], 
[S/H], [Si/H], [Fe/H]+0.4 (where the offset accounts for typical 
dust depletions; Prochaska & Wolfe 2002). Since the correction to 
[Fe/H] is rather crude we give priority to, for example, a [Zn/H] 
limit over a [Fe/H] detection. The exception to these priorities is 
AO 0235+ 164 for which we use an X-ray metallicity, as there is 
evidence that this is a very dusty absorber with depletion consid- 
erably larger than the typical value of 0.4 that we otherwise adopt 
(e.g. Junkarrinen et al 2004; Kulkarni et al. 2007). The lower left 
part of the diagram with JV(Hi)< 2 x 10 20 cm" 2 and T,< 600 K 
is poorly populated by the Galactic and sub-DLA sightlines alike. 
This supports the suggestion by Kanekar et al. (2011) that below 
this column density threshold, the neutral medium is dominated by 
warm gas. At higher column densities, the DLAs show much more 
scatter in T, than the Galactic data, but with a clear dependence 
on metallicity. The DLAs with metallicities below around 1/10 of 
the solar value (yellow, green and blue points) have high spin tem- 
peratures for their Af(Hi), relative to the Milky Way points. The 
more metal-rich DLAs (orange and red points) have either compa- 
rable, or low spin temperatures for their A'(Hi). Such a metallic- 
ity dependence may be expected from cooling arguments, whereby 
metals aid the cooling of the neutral medium (Wolfire et al. 1995; 
Kanekar & Chengalur 2001), in a similar way that metals have been 
suggested to aid the transformation from the neutral to molecular 
regimes (Schaye 2001; Krumholz et al. 2009). 

In Figure [8] we make a crude conversion of the spin temper- 
ature into a CNM fraction (see also Kanekar & Chengalur 2003). 
For this calculation, in which we adopt a two phase medium, we 
must assign spin temperatures to the cold and warm neutral me- 
dia (WNM). For the CNM, the spin temperature is the same as 
the kinetic temperature (e.g. Liszt 2001), and we adopt a value of 
9OK0. However, in the WNM, the low densities are insufficient to 
thermalise the Hi 21 cm transition (e.g. Liszt 2001). Therefore, al- 
though the kinetic temperature of the WNM is expected to be 5000- 
8000 K, the spin temperature, which will depend on the physical 
conditions within the cloud, is expected to be considerably lower, 
1000-4000 K (Liszt 2001). Although the CNM fractions that are 
calculated will depend on the temperatures used, the qualitative dis- 
tribution of points remains the same. We adopt a T, for the WNM 
of 3000 K, although we note that this will lead to an (unphysical) 
negative CNM fraction if the measured T s > 3000 K. Figure[8]shows 
that, as suggested by Kanekar et al. (201 1), the low A^H i) Galactic 
sightlines are dominated by the WNM. As the A'(Hi) increases, 
a higher CNM fraction is present. The sub-DLA in our sample 
also has a relatively low CNM fraction, consistent with the upper 
end of the distribution seen in Galactic (sub-DLA) sightlines. In 
the DLA regime, the low metallicity absorbers have relatively low 
CNM fractions, compared with the Milky Way. The DLAs with 
high CNM fractions tend to be relatively metal-rich. 



4 Although CNM temperatures of 100-200 K are often assumed (e.g. 
Kanekar & Chegalur 2003), the lowest T, in our sample is 90K. So as not 
to exceed CNM fractions of unity, we adopt 90 K as the CNM temperature. 
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Figure 6. Abundances taken from the literature for DLAs and sub-DLAs (red and blue open circles respectively). For clarity, limits are not shown for the 
literature sample of DLAs and sub-DLAs for which there are no Hi 21 cm observations. Hi 21 cm absorbers with measured covering factors are shown with 
stars, again coloured coded accorded to their classification as DLAs (red) or sub-DLAs (blue). Open stars show H 1 21 cm detections and limits with T, > 500 K 
and filled stars are for H i 21 cm absorbers with T, < 500 K. 
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Figure 7. Covering factor-corrected spin temperatures versus A'(Hi) for 
DLAs, colour-coded by metallicity ([Zn,S,Si,Fe+0.4/H]). PDLAs are ex- 
cluded from this Figure. For comparison, Galactic sightlines taken from 
Kanekar et al. (2011) are shown with smaller black points. The horizon- 
tal dashed line shows a spin temperature of T s = 600 K. The vertical 
dashed line indicates the division between sub-DLAs and DLAs (at log 
iV(Hi) = 20.3). 




Log N(HI) 

Figure 8. The fraction of gas in the CNM versus N(H i) for DLAs, colour- 
coded by metallicity ([Zn,S,Si,Fe+0.4/H]). PDLAs are excluded from this 
Figure. For comparison, Galactic sightlines taken from Kanekar et al. 
(2011) are shown with smaller black points. The vertical dashed line in- 
dicates the division between sub-DLAs and DLAs (at log iV(Hi) = 20.3). 

3.4 The spin temperature-metallicity anti-correlation 

In the previous sub-section we showed that sub-DLAs tend to have 
high T s (and therefore high fractions of warm gas), regardless 
of their often elevated metallicities. DLAs may exhibit significant 
fractions of cold gas, but this is likely to be metallicity dependent. 

We now turn to the explicit dependence of spin temperature 
on metallicity. Kanekar & Chengalur (2001) predicted that an anti- 
correlation between T, and metallicity would naturally arise if the 
paucity of cold neutral gas in DLAs was due to fewer cooling routes 
due to a lack of metals in the absorbers. The largest sample to have 
been previously used to address this issue was a sample of 26 ab- 
sorbers, of which 10 systems had both T, and metallicity measure- 
ments and 20 systems had estimates of the low-frequency covering 
factors (Kanekar et al. 2009c). The number of detections in our 
sample is significantly increased to 17, even with the additional 
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Figure 9. Spin temperature as a function of metallicity 
([Zn,S,Si,Fe+0.4/H]) for all DLAs (red) and sub-DLAs (blue) with 
measured covering factors. Absorbers from our new sample (this paper) 
are circled. The open star represents AO 0235 + 164 for which we adopt the 
X-ray metallicity from Junkkarinen et al. (2004). The upper panel includes 
metallicity and/or spin temperatures limits. The lower panel includes only 
absorbers for which both H i 2 1 cm and metal absorption is detected. The 
dotted line shows the best fit for the detections taken from Kanekar et al. 
(2009), whereas the solid line shows the best fit for the full data sample 
presented here 
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requirement that all absorbers must have low frequency covering 
factor measurements to be included in our sample. The new data 
presented here include DLAs with very high and low spin temper- 
atures, which are particularly useful for anchoring any relationship 
with metallicity (Figure [9). We note that the lowest metallicity ab- 
sorber in the lower panel of Figure [9] (which has the highest spin 
temperature amongst the systems with both Hi 21 cm and metal 
line detections) is Q0952+179. Due to our prioritization of metal- 
licities, this DLA appears in the upper panel as an upper limit from 
[Zn/H], but its measured value of [Fe/H] is used in the lower panel. 

For the full sample of 26 detections and limits, a non- 
parametric generalised Kendall-tau rank correlation test (the BHK 
statistic in the ASURV package; Brown et al. 1974; Isobe et 
al. 1986) detects an anti-correlation between metallicity and spin 
temperature at 3.4cr significance (Figure [9] upper panel), treating 
the metallicity as the independent variable. If the DLA towards 
AO 0235+ 164 is dropped from the sample (in contrast with the rest 
of the sample, its metallicity has been obtained from a X-ray mea- 
surement), the significance of the anti-correlation is 3.1cr. Reduc- 
ing the sample to only the 17 systems with measurements of both 
T, and metallicity (i.e. excluding limits; see Figure[9] lower panel) 
yields an anti-correlation with 3. Oct significance. Again, dropping 
AO 0235+ 164 from this detection sample gives an anti-correlation 
with a 2.8<r significance (16 systems). 

There is only one sub-DLA included in Figure [9] the za> B = 
1.6724 absorber towards Q0237-233, whose spin temperature we 
present for the first time in this work. Despite the low ,/V(H i), the 
spin temperature is relatively low in this absorber: T ,= 390± 127 K. 
The moderately high metallicity (Z ~ 1/4Z Q ) places this sub- 
DLA in good agreement with the DLA data on Figure [9] Drop- 
ping the one sub-DLA from the sample yields a T,-metallicity anti- 
correlation that is significant at 3.2<r and 2.7cr for limits+detections 
and detections only, respectively. 

We also used a linear regression analysis to obtain the best-fit 
relation between spin temperature (log[T s ]) and metallicity [Z/H], 
applying this to the 17 DLAs with measurements of (rather than 
limits on) both T s and [Z/H] (see Figure [9] lower panel). We used 
the BCES(Y|X) estimator (Akritas & Bershady 1996) for this pur- 
pose, again treating [Z/H] as the independent variable X. Note that 
this method takes into account measurement errors on both vari- 
ables, as well as correlations between these errors. This is relevant 
here as both T s and [Z/H] are derived from the Hi column den- 
sity, and the errors on the two quantities are hence correlated. The 
regression analysis yielded log[T s ] = (-0.87 ± 0.14) x [Z/H] + 
(1.93 ± 0.19), consistent with the results of Kanekar et al. (2009c). 
The fit is shown as a solid line in the lower panel of Figure [9] 
If, instead of the (Y|X) technique, we adopt the BCES bisector 
method for independent variables, the fit is consistent to within ltr: 
log[T s ] = (-1.11 + 0.16) x [Z/H] + (1.70 + 0.19). 

Kanekar et al. (2009c) did not find a significant anti- 
correlation between dust depletion and T,. In Figure [10] we confirm 
this by plotting T s versus dust depletion [X/Fe] (where X=Zn, S, 
Si, or the X-ray value for AO 0235 + 164) for our expanded sample. 
Applying the generalised Kendall tau test to our sample of 22 sys- 
tems with either measurements of or limits on T s and dust depletion 
(FigureQI)] upper panel) we find a 1 .9<x anti-correlation between T s 
and [X/Fe], considerably weaker than the anti-correlation between 
T s and metallicit}0. Note that only a handful of DLAs (11 systems) 
have measurements of both dust depletion and T s ; this sub-sample 

5 Adding a further two absorbers to the depletion sample which do not 
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Figure 10. Spin temperature as a function of dust depletion as measured 
by [X/Fe]=[Zn,S,Si/Fe] for absorbers with measured covering factors. The 
upper panel shows detections and limits in H i 21 cm and metals. The lower 
panel shows only detections in Hi 21 cm and metals. Sub-DLAs are shown 
in blue, DLAs are shown in red and absorbers from our new sample (this 
paper) are circled. PDLAs are excluded from this Figure. 

shows no evidence of an anti-correlation between [X/Fe] and T s , 
with a significance of ~ 0.4cr. The lack of a significant depen- 
dence of T, on depletion demonstrates that the T,-metallicity anti- 
correlation is not driven by an underlying dependence on dust. We 
also note that the DLAs with Hi 21 cm absorption show a range of 

have Fe abundances, but have detections or limits in Cr n increases the sig- 
nificance negligibly to 2.0 cr. 
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Figure 11. Metallicity, as traced by [Zn/H] as a function of N(H i) for a liter- 
ature sample of optically selected DLAs (black open circles) and Hi 21 cm 
absorbers (red stars). Open stars show H i 21 cm detections and limits with 
T,> 500 K and filled stars are for H i 21 cm absorbers with T s < 500 K. 



depletions that is well-matched to the range in the full DLA sample, 
where the mean value of [Zn/Fe] — 1-0.4 (e.g. Prochaska & Wolfe 
2002). 



3.5 H i 21 cm absorption as a tool for selecting DLAs. 

The absorbers in the sample selected here were identified first as 
Mgn systems, then as Hi 21 cm absorbers, with the A'(Hi) deter- 
mination as the final step. Excluding B0801+303, which was post 
facto found to be an extended radio source, 4/5 of the Mg n-selected 
H 1 21 cm absorbers in our sample are found to be DLAs. This com- 
pares with a 35% DLA identification rate for absorbers selected 
purely on Mg n equivalent width (Rao et al. 2006) and ~ 60% when 
velocity widths are also taken into account (Ellison 2006; Ellison, 
Murphy & Dessauges-Zavadsky 2009). The one sub-DLA in our 
main sample is towards B0237-233, whose Hi 21 cm optical depth 
is the lowest in our sample by a factor of ~ 3. 

The connection between H i 21 cm absorbers and DLAs is not 
particularly surprising. Several works (e.g. Kanekar & Chengalur 
2003; Curran et al. 2010) have noted a correlation between the in- 
tegrated H i 21 cm optical depth and N(H i). Further, Kanekar et al. 
(2011) show that most Galactic sightlines with N(Hi)< 2 x 10 20 
cnT 2 have relatively high spin temperatures, and are hence domi- 
nated by the WNM. Only above this approximate threshold can the 
neutral gas apparently cool sufficiently to form a significant fraction 
of the cool phase, then becoming readily detectable in Hi 21 cm ab- 
sorption. 

Hi 21 cm absorbers may therefore represent an interesting al- 
ternative identification of DLAs at low redshifts. The obvious con- 
cern is how representative such a selection would be. Given the 
results discussed in this paper, it may be expected that Hi 21 cm 
absorption is most readily detected at large N(H i) and high metal- 
licity. Indeed, all four DLAs of the present sample have high N(H i), 
in excess of 10 21 cm 2 . This is unsurprising, since for a given T s 
the Hi 21 cm optical depth will be highest for the highest A'(Hi) 
absorbers. Since DLAs with log A'(Hi) > 21 are relatively rare, 
H 1 21 cm selection offers a good opportunity to identify the highest 
column density absorbers for targeted studies. For example, many 
of the 'metal-strong' DLAs (Herbert-Fort et al. 2006; Kaplan et al. 



2010), which are useful for the identification of rare metal species, 
have such high A'(Hi). 

Despite the bias towards high A'(Hi) in an Hi 21 cm selected 
sample of DLAs, the metallicities and dust depletion in Hi 21 cm 
absorbers span the observed range in the full sample of DLAs and 
sub-DLAs (Fi gures |6land| 1 0| . In FigureQTJ we also show the metal- 
licities (as measured by [Zn/H]) as a function of N(K i). It has been 
argued that the anti-correlation between [Zn/H] and iV(Hi) is in- 
consistent with a dust bias, given the typical reddening in DLA 
samples (Elllison, Hall & Lira 2005). The paucity of high N(H i) 
absorbers with high metallicities can instead be understood as a 
metallicity-dependent transition from the atomic phase to a phase 
dominated by molecular gas (Schaye 2001; Krumholz et al. 2009). 
The Hi 21 cm absorbers have metallicities consistent with those of 
the full sample of DLAs at a given N(H i). Combined with Figures 
[6l and [Tol this indicates that an Hi 21 cm selected sample is not 
strongly biased for studies of elemental abundances. 



4 CONCLUSIONS 

The study of Hi 21 cm absorbers at cosmological distances is now 
coming of age in terms of sample sizes, redshift coverage and avail- 
ability of covering factors. This paper focuses on our efforts to fill 
in the spin temperature redshift desert (0.6 < z < 1.7) and investi- 
gate the role of metallicity in modulating spin temperature trends. 
In a companion paper (Kanekar et al. in prep.) much of the support- 
ing high redshift data used in our analysis is presented, as well as a 
analysis of the redshift evolution of the spin temperature. The main 
results of the current paper are as follows. 

(i) For a sample of six Mg n-selected Hi 21 cm absorbers in the 
redshift range 0.6 < z a b s < 1.7 we have determined A'(Hi) (mostly 
from HST spectra), metal column densities from optical echelle 
spectra and covering factors (/) from VLBA images. 

(ii) One of the absorbers (z abs = 0.60190 towards J1431+3952) 
is a new H i 21 cm detection made with the GBT Its spin tempera- 
ture is the lowest yet reported: T,=90±23 K. 

(iii) The metallicities and dust depletions of H 1 21 cm absorbers 
span the range that is typical for the full DLA sample in the 
same redshift range (or at fixed A'(Hi)). As is the case for the full 
DLA, Hi 21 cm sub-DLAs have slightly higher metallicities than 
Hi 21 cm DLAs, at a fixed redshift. 

(iv) We confirm the presence of an anti-correlation between spin 
temperature and metallicity, using an absorber sample with mea- 
surements of the covering factor. This is detected at 3.4cr signifi- 
cance for the sample containing both measurements and limits of 
T, and metallicity (26 systems) and at 3.0cr significance for a sub- 
sample (17 systems) with measurements of both T, and metallicity. 

(v) Although a wide range of spin temperatures are found at 
fixed N(H i), the CNM fraction clearly depends on metallicity. Ab- 
sorbers with low metallicity (for their N(H i)) tend to have high spin 
temperatures and low CNM fractions. Only the moderate-to-high 
metallicity (above ~ 1/10 solar) DLAs tend to show CNM fractions 
comparable to, or in excess of, Galactic sightlines, at fixed neutral 
hydrogen column density. 
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APPENDIX A: METAL LINES OBSERVED IN OPTICAL 
ECHELLE SPECTRA 

Examples of the echelle data obtained either from the UVES 
archive or obtained by us from HIRES are presented in Figures 
lAT1tolA5l 
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Figure Al. Metal line transitions in the DLA towards B0105-008. The velocity scale is relative to z = 1.37104. 
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Figure A2. Metal line transitions in the DLA towards Q0237-233. The velocity scale is relative to z = 1.67235. 
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Figure A3. Metal line transitions in the DLA towards J1431+3952. The velocity scale is relative to z = 0.6018. 
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Figure A4. Metal line transitions in the DLA towards J1623+0718. The velocity scale is relative to z = 1.3357. 
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Figure A5. Metal line transitions in the DLA towards B2355-106. The velocity scale is relative to z = 1.17230. 
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